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Abstract. We present results of an observation of the anomalous X-ray pulsar 4U 0142+61 with the EPIC cameras on XMM- 
Newton performed on 2003 January 24. The pulse phase averaged spectrum can be best described by the sum of a black 
body with a temperature of kT^^ = 0.395(5) keV and a power law with photon index F = 3.62(5). The unabsorbed 2-10 keV 
flux is ~7.2 • 10""ergs"' cm"'. No evidence for additional spectral features such as cyclotron lines is present. These results 
are consistent with those from an earlier Chandra observation in 2001 May. Phase resolved spectroscopy over the 8.6882(2) s 
period (MJD 52663.93) shows clear variations with pulse phase of F, while A:rBB shows a small variation of ~ 12%. We confirm 
earlier conclusions by Qzel (.2001) that the emission from AXPs is more likely to originate from neutron stars with one single 
and magnetically active region on the neutron star The significantly diflrerent behavior of the soft and hard spectral components 
with pulse phase, however, cannot be fully reconciled with the present magnetar emission models. 
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1. Introduction 

4U 0142+61 is a member of the small group of anomalous X- 
Ray pulsars (AXPs) which are characterized by spin periods in 
the range of 8- 12 sec, associated with a relatively stable long 
term spin down ('Gavriil & Kaspi"2002'), very soft X-ray spectra 
^ well modeled by blackbody plus power law components with 
^ r ~ 3, a X-ray luminosity of ^lO^'^-lO^^erg s~\ l ow in com- 



pariso n to High Mass X-Ray Binaries (HMXB) jlsrael et alJ 
I2OO2I) . and very faint optical counterparts. It is generally as- 
sumed that these objects are isolated neutron stars as no binary 
companion for any of these objects has been found up to now. 

The X-ray luminosity of AXPs exceeds the energy avail- 
able from the spin down of the neutron star, so some additional 

1 • ^1 1 

energy source is needed ( van Paradiis et al. 1995). There are 

mainly two classes of models trying to explain the source of 
the missing energy. One class proposes that the AXPs are pow- 
ered by accretion, either froi n debris of a di srupted HMXB after 
a common envelope phase Jvan Paradiis et al. 1995 ) or from a 
disk formed by the fallback materi al from a supernova explo- 
sion ( Chatteriee & Hernauist '2000^. The other class of models 
assumes a very strong magnetic field which can explain the 
spin down as magnetic dipole radiation while the X-ray lu- 
minosity is supplied by either the decay of the magnetic field 

1 Thompson & Dutican 1995) or the cooling of the neutron star 
Hevl & HernauisJl997l) 

The brightest of the small sample of anomalous X-ray 
pulsars is 4U 0142+61, a source detected by Uhuru, which 



had its nature established during an EXOSAT campaign in 
which a 8.7 s periodic variation was found (Israel et al. 193)- 
Subsequent X-ray observations manifested features usually at- 
tributed to AXPs such as a spectrum which is best fit by a 
0.386(5) keV bl ack body p l us a p ower law with photon index 
r = 3.67(9). IWhite et al.l ( Il996l) . Long term period changes 
were discovered in subsequent ASCA and RXTE observa- 
tions (Pauletal. 2000; Gavriil & Kaspi 2002), with the lat- 
ter revealing a period of P = 8.68832877(3) s and a period 
change of P = 2.02(25) x 10"'^ s s"'. Finally, an optical coun- 
terpart was detected which s howed strong p eriodic variation 
dKern & Martin' 2002. see also lHuUeman et alJ2004 . A recent 
observation of 4U 0142+61 by Chandra confirmed the associ- 
ation of this optical counterpart with the X-ray source and also 
resulted in a further improvement of the X-ray spectral results 
JPatel et al.l2003Hjuett et alJ2002l) . 

In this Paper we report an observation of 4U 0142+61 with 
XMM-Newton. In Sect.|2lwe describe the observation and data 
reduction. Lightcurve analysis and pulse phase averaged spec- 
tral analysis, presented in Sect. |3l and Sect. |3 confirm ear- 
lier Chandra and RXTE results. Phase resolved spectroscopy 
shows a variation of two different spectral components which 
are presented in Sect.|5]and discussed in Sect.|6l 

2. Observations witli XMM-Newton 

The observation of 4U 0142+61 with XMM-Newton was per- 
formed on 2003 January 24 (XMM-Newton Revolution 573). 
This observation has an exposure time of ~4.2ksec with the 
EPIC-pn camera in small window mode and ~5.8ksec with 
both EPIC-MOS cameras in fast timing mode. The EPIC-pn 
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detector is a back-illuminated charge coupled device (CCD) us- 
ing 6 4 X 64 pixels in sma ll window mode, covering a ~4' x 4' 
area jStriideret al1l200ll) . The two MOS cameras are front- 
illuminated CCDs detectors jTurner etalJ200ll) . In their timing 
mode, spatial information is preserved only along the x-axis of 
the detector, but events are registered with a higher time reso- 
lution than in the EPIC-MOS imaging modes. 

Standard data processing was performed with the XMM- 
Newton Scientific Analysis Software, version 5.4.1. For the 
EPIC-pn the source region was taken to be a circle of radius 
~ 100" around the point source. The background was measured 
in two rectangular regions arranged at the corners of the read- 
out window. For the MOS cameras, events in the raw x interval 
280 . . . 330 were chosen for source events, background events 
were taken from the intervals 260 . . .279 and 331 .. . 345. No 
background flares were present during this observation. We 
checked the EPIC-pn data for evidence for pileup and found 
no indication for pattern pileup or pileup in diagonal pixel pat- 
terns at energies above 0.3keV. Fo r the MOS 1/2 da ta pattern 
pileup is negligible in timing mode llEhle et al.l2003l) . 

We note that a second earlier XMM-Newton observation 
of 4U 0142+61 exists, which was performed on 2002 February 
13 (revolution 399). For this 3.4ksec long observation, only 
data from the EPIC-pn cameras are available (the MOS cam- 
eras were switched off) and the background was a factor of 10 
higher than in the Rev. 573 observation, such that we decided 
to not use these data here. 

The position of 4U 0142+61 was estimated by fitting a 
Gaussian to the point-spread function of the EPIC-pn data 
in both observations. The position obtained is Q'j2()oo.o = 
01^ 46" 22'.'3, (5j2()oo.o = 61°45'02'.'5, which is consistent within 
the XMM 2'.'2 error circle with the more accurate Chandra po- 
sition of a = l''46'"22?42, 6 ^ 6r45'02'.'8 (positional un- 
certainty ~ y.^ JPatel et alJ2003h . 

3. Timing analysis 

For all instruments the lightcurve was s earched for period s 
in the 8.6 s domain using epoch folding llLeahv etalJ ll983V 
Combining the results of the period search from all instru- 
ments of Rev. 573 (Epoch: MJD 52663.93) we found a period 
of P ^ 8.6882(2) s. This period agrees with the RXTE data 
of lGavriil & Kaspi (2002) which predict P = 8.68849(1) s, as- 
suming a constant P = 2.02(25) x 10"'^ s s"' . 

We estimate the uncertainty of the XMM-Newton period 
using a Monte-Carlo simulation: A synthetic lightcurve was 
constructed by taking the time stamps as in the measured 
light curve and associate them with the flux of the folded 
pulse profile at the corresponding phase (for the best period 
and ephemeris found from the observed data). From this tem- 
plate 1000 simulated lightcurves were constructed assuming 
the mean pulse profile and then randomizing the flux values 
using the appropriate Poisson distribution. For all simulated 
lightcurves the best periods were determined. The standard de- 
viation of the distribution of those periods was taken as the 
uncertainty of the period found in the observed data. 

For further timing and profile analysis we used a period of 
8.68823 s. To compute pulse profiles all single events from all 
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Fig. 1. Pulse profile of 4U 0142+61 in the energy range 0.3- 
lOkeV. 
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Fig. 2. Profile of the combined lightcurves in the energy ranges 
0.3-1. 3keV, 1.3-2.5keV, and 2.5-lOkeV. The peak-to-peak 
pulse fractions are: 8.6% + 1.3% (0.3-1. 3 keV), 9.1% + 0.9% 
(1.3-2.5 keV), and 9.6% +2.0% (2.5-lOkeV). The dash dotted 
line is the result of the fit of a gaussian distribution to the peak 
at phase ~ 0.9. Here all data from EPIC-pn and MOS-1/2 from 
Rev. 573 were used. 



three EPIC detectors are added. Multiple events were not taken 
into account. The pulse profile in the energy range 0.3-lOkeV 
is shown in Fig. [2 Its peak-to-peak pulsed fraction, defined as 

(i^max - /^mm)/(^'max + i^mm), is 7.7% + 0.9%. 

In Fig. |2] we display the pulse profile and give the peak- 
to-peak pulse fractions for several energy bands. As was first 
noted bv lWhite et akld 19961) the pulse is double peaked at low 
energies. At higher energies the main pulse remains as a broad 
single peak centered at phase ~0.9, while the secondary peak 
weakens and changes into a double peaked feature. 

In our data we observe a slight phase shift of the main peak. 
Fitting a Gaussian distribution to this peak yields a centered 
phase value of 0.878(9) (0.3-1. 3keV), 0.891(5) (1.3-2.5 keV) 
and 0.92(1) (2.5-lOkeV). The peak-to-peak pulsed fraction is 
not significantly energy dependent, with measured pulsed frac- 
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Fig. 3. Spectrum and residuals for the blackbody plus power 
law model. 

tions of 8.6% + 1 .3% (0.3-1 .3 keV), 9. 1 % +0.9% ( 1 .3-2.5 keV), 
and 9.6% ± 2.0% (2.5-lOkeV). 

4. Spectra 

For spectral analysis we used the single event data from 
the EPIC-pn camera and the MOS cameras of Rev. 573. 
Simultaneous fits were made with XSPEC version 11. 2. Obi in 
the energy ranges of 0.5-9.5 keV (pn) and 0.5-8 keV (MOS), 
using a multiplicative constant to take into account the flux cal- 
ibration uncertainty of the MOS cameras with respect to the 
EPIC-pn. To allow analysis the channels were rebinned to 
contain at least 50 counts per bin. 

Several models were used to fit the phase averaged spectra: 
single component models (blackbody, bremsstrahlung, power 
law) and combinations of those, including the classical combi- 
nation of a black body component plus a power law. All models 
include absorption due to the interstellar medium with hydro- 
gen equivalent column density A^h (using the XSPEC model 
phabs). Fit results of the different models are shown in Tabled 

Single component models gave unacceptable results. The 
best fit is obtained for the black body plus bremsstrahlung 
model. The blackbody plus power law model (Fig.^} gave al- 
most the same x'' but had a smoother residual at higher ener- 
gies, and we will use this model in our further discussion. The 
residuals of this model show no indication for any line in the 
energy band covered by our observation. Using this model the 
unabsorbed 2-lOkeV flux is 7.02 x 10"'' erg s"' cm"^. No fur- 
ther spectral features such as cyclotron lines were found in the 
spectral range of 0.5-10 keV. 

Extending the fit to energies below 0.5 keV for the EPIC-pn 
data alone and refitting yields a deviation of ~30% between the 
data and the model in 0.3-0.5 keV band. This deviation i s larger 
than allowed by the uncertainties of the calibration jKirschl 
I2OO3I) . Phase resolved spectroscopy shows this excess of the 
model in all phases (see Sect. |5}, i.e., the excess flux is not 
caused by a possible continuum variability with pulse phase. 
The feature also is present for different absorption models and 
abundances and is also present in the double event spectrum of 
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Fig. 4. Results of phase resolved spectroscopy for a blackbody 
plus power law model. Shown are a total EPIC count rate, b 
blackbody temperature, c photon index, d blackbody flux, and 
e flux of the power law. The 0.5-lOkeV fluxes are unabsorbed 
fluxes and are in units of 10 ' ' erg cm"^ s ' . 

the EPIC-pn camera. We assume that the residual either comes 
from a higher galactic absorption than the absorption model or 
from unresolved calibration issues in this energy band. For the 
remainder of this paper, therefore, we will restrict ourselves to 
the energy band above 0.5 keV. 

5. Phase resolved spectroscopy 

To search for spectral variations with pulse phase we split the 
data into ten phase bins, rebin the spectra to contain at least 25 
counts per bin, and model the EPIC data of each phase bin with 
a black body plus power law model. Fig.|4t shows the resulting 
pulse profile. 

In a first step, all parameters were allowed to freely vary 
with pulse phase. We found that the column density A^h did 
not show significant variation. To test this constancy we fit the 
subset of phase bins with the highest flux, constraining A^h to be 
the same for these bins. In a second iteration, A^h was allowed 
to vary with pulse phase. The F-test shows a probabihty of 
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Table 1. Phase averaged spectral fits. Best fit parameters for an absorbed power law (PL) and a bremsstrahlung (BREMS) con- 
tinuum, with and without an additional blackbody (BB) component. The parameters are the equivalent hydrogen column density. 
Nil, the power-law photon index, F, and its normalization at 1 keV in units of 10"' Photons keV^'cm^^s"', the temperature of the 
bremsstrahlung spectrum, ^rbiems. and its 2-lOkeV component flux contribution, Fbiems^ in 10""ergs"' cm"^ for the EPIC-pn 
camera, the temperature of the black body, kT},^, and the radius of the black body, /?|,b, assuming a distance of lOkpc. MOS 1/2 
spectra are scaled with constants S'mosi and 5mos2 with respect to the EPIC-pn. All error bars are at the 90% level for one 
interesting parameter and expressed in parentheses in units of the last digit shown. 
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7 1 % that A^H does not vary. Therefore we considered A^h to be 
a constant parameter with a value of Ah - 0.91 x 10^^ cm"^ in 
our subsequent analysis. 

The results of our spectral analysis are shown in Figs.|3) 
and c as a function of pulse phase. The blackbody temperature 
is virtually constant, while the photon index shows a clear pulse 
phase dependence. Using the unabsorbed spectrum, we further- 
more derive the contribution of the blackbody and the power 
law component to the total 0.5-lOkeV flux as a function of 
pulse phase. The result is shown in Figs.|3l and e. Uncertainties 
are computed using a Monte-Carlo approach randomizing the 
input data sets. We compare these results with AXP emission 
models in the following section. 

6. Discussion 

In this paper we have shown that the comparatively short obser- 
vation of 4U 0142+61 with XMM-Newton confirms the spec- 
tral and timing parameters from earlier observations. The best 
fit models are the standard black body plus power law and the 
black body plus bremsstrahlung model. The period agrees well 
with the high precision value derived from the long term timing 
campaign with RXTE iGavriil & Kaspi' 2002'). 

The flux of 7.02 x IO~"erg sec ' cm"^ is slig htly lower 
than th e flux reported in the ASCA observations of IPauI et"al] 
ranging from 9 ... 16 x 10 "erg sec ' cm ^. We con- 
sider this discrepancy as negligible because the authors of the 
ASCA observation note the possibility of a flux overestimation 
due to limited spectral resolution. 

Pulse phase resolved spectroscopy allowed us to disentan- 
gle the different phase dependencies of the relevant parame- 
ters (Fig.E): The blackbody temperature is constant with pulse 
phase, although taken at face value there is a ~12% variation 
with phase and a slightly lower temperature at phase ~0.3. The 
flux of both, black body and power law, follows the double 
peaked pulse profile (Fig.|3t), although but for the powerlaw 
flux variation the peak at phase ~ 0.5 is suppressed. This is 
reasonable if one takes into account the hardening of the power 
law at phase ~ 0.5. 

As has been pointed out, e.g., bv lOzell ( l200lb and lOzell 
il2002h . observed pulse profile variations in AXPs can be used 
to test emission models for magnetars. These models predict 



a fan beam plus pencil beam emission pattern and a resulting 
pulse profile variation that is strongly dependent on the charac- 
teristics of the emitting plasma. The observed pulsed fraction 
depends on the angle between the magnetic field axis and the 
rotation axis of the neutron star as well as on the angle between 
th e observer's line of sight and the rotational axis. As shown 
bv lOzel et al] (I2OOI ), the high pulsed fractions of most AXPs 
argue for the radiation to come from a single emitting region on 
the neutron star, e.g., cause d by a region of enhanced magnetic 
activity on the neutron star ( I6zelll2003) . although the parame- 
ters of 4U 014 2+61 are sti ll barely compatible with two emis- 
sion regions ( Ozel et al.l2001il) . We can confirm this conclusion 
of a single emission region from the XMM-Newton data, which 
shows only a small change of the pulsed fraction with energy, 
which is in contrast to the strong energy dependence expected 
from models with two antipodal emission regions. 

Using phase resolved X-ray spectroscopy, it is in principle 
possible to further constrain these emission models. In terms 
of the empirical spectral model applied here, magnetar models 
predict that the blackbody and powerlaw components should 
be coupled through resonant Compton heating of photons scat - 
tered at cyclotron resonances (see also ' Thompson et a Dl2002h . 
This coupling results in pulse profiles which are predicted 
to have essentially the same shape over the 1-lOkeV energy 
range. In our observations, this predicted behaviour cannot be 
fully confirmed: Both, the energy resolved profiles as well as 
the decomposition of the pulse into a blackbody and a pow- 
erlaw component show only little coiTelation. For example, 
Pearson's correlation coefficient, r, of the component fluxes of 
the black body and the powerlaw is only r - 0.561. On the 
other hand, however, the similar variation of both model com- 
ponents suggest that there is a common flux contribution super- 
imposed with an spectral variation of the powerlaw. This would 
explain the decrease of the pulse at phase ~ 0.5 with higher en- 
ergies. 

In summary, in this paper we were able to show that phase 
resolved spectroscopy adds power to constrain existing models 
and allows in principle to test different hypotheses describing 
the physical geometry of the system. Further and longer obser- 
vations and further modeling are required, however, until these 
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tests will yield a conclusive and selfconsistent answer on the 
emission process at work. 
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